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POTENTIAL FLOW ABOUT ELONGATED BODIES OF REVOLUTION 

By Carl Kaplan 



SUMMARY 

This report presents a method of solving the problem 
of axial and transverse potential flows around arbitrary 
elongated bodies of revolution. The solutions of La- 
place's equation for the velocity potentials of the axial 
and transverse flows, the system of coordinates being an 
elliptic one in a meridian plane, are known to be of the 
following form: 

CO 

n=i 

03 

<j> = ^^AmQn 1 (X) P n l 0*) cos 0 (transverse flow) 

n=i 

If a power-series development of X in m is assumed 
as the equation of the meridian profile in elliptic co- 
ordinates, the boundary conditions of the two types of 
flow yield linear equations for the determination of the 
coefficients A n and A n \. It is further shown that a 
knowledge of these coefficients leads directly to the sink- 
source and doublet distributions corresponding to the 
axial and transverse flows, respectively. 

The theory is applied to a body of revolution obtained 
from a symmetrical Joukowsky profile, a shape re- 
sembling an airship hull. The pressure distribution 
and the transverse-force distribution are calculated and 
serve as examples of the procedure to be followed in the 
case of an actual airship. A section on the determina- 
tion of inertia coefficients is also included in which 
the validity of some earlier work is questioned. 

INTRODUCTION 

There are two methods of handling the problem of 
potential flow about a body of revolution. One, the 
indirect method first used by Taylor (reference I) 
and by G. Fuhrmann (reference 2) who computed the 
pressure distribution by the method of sources and 
sinks suggested by Rankine. Fuhrmann assumed 
certain sink-source distributions and calculated the 
pressure distribution for the streamline body resulting 
from the assumed sink-source system. He also con- 
structed models of the calculated shapes and measured 



the pressure distributions over them when placed in a 
wind tunnel. 

The other method, developed by von Karman 
(reference 3), considered the direct problem; i. e., the 
calculation of the pressure distribution over a given 
streamline shape. He approximated the requisite 
sink-source distribution by a computed continuous 
system of sinks and sources arranged in stepwise con- 
stant intensity. The various strengths were deter- 
mined from the condition that the airship hull is a 
streamline surface in the parallel flow and the flow 
induced by the sinks and sources. By satisfying this 
condition at an arbitrary number of points equal to the 
number of unknown sink and source segments, von 
Karman obtained a system of linear equations for the 
determination of the unknown strengths of the sink- 
source distribution. He also treated the case of 
transverse flow (references 3 and 4) by the distribution 
of doublets along the axis of symmetry of the body of 
revolution and calculated the strengths of the various 
doublet segments in a manner similar to that used for 
the sink-source intensities. 

The present paper is an attempt to treat the direct 
problem according to the methods of the potential 
theory. Thus, Laplace's equation for the velocity 
potential is set up in a system of elhptic-cylindrical 
coordinates X, p., 8 and solved in conjunction with the 
appropriate boundary conditions for axial and trans- 
verse flows. It is then assumed that a power-series 
development of X in p represents the meridian profile 
of the elongated body of revolution. The boundary 
conditions for the two types of flow may then be ex- 
pressed in the form of power series in n valid for the 
entire range of /i. This method leads to two sets of 
linear equations, each set infinite in number of equa- 
tions and each equation containing an infinite number 
of unknown coefficients which serve to determine the 
velocity potentials for the axial and transverse flows. 
Incidental to the major task of deterniining these 
coefficients, the sink-source and doublet distributions 
corresponding to the axial and transverse flows are also 
determined. Thus the results of this method are essen- 
tially the same as those obtained by the method of 
von Karman but are obtained in a more rigorous and 
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direct m ann er. In von Karman's method, approxima- 
tions are made prior to the analysis; whereas, in the 
method presented in this paper, approximations are 
made after the analysis has been carried through in a 
rigorous maimer. 

FUNDAMENTAL EQUATIONS 

The fluid motion is assumed to he steady and irrota- 
tional. There then exists a velocity potential <f>, 
which is, in general, a function of the rectangular 
Cartesian coordinates (x, y, z). In cases of rotational 
symmetry, however, it is appropriate to introduce the 
cylindrical coordinates (z, p, 8) where z denotes the 
distance along the axis of symmetry, p(=V^+1/ 2 ) tne 
perpendicular distance from this axis, and 8 the angle 
between the (z, p) and (z, a:) planes. (See fig. 1.) 




Figure 1. 



Furthermore, since only elongated surfaces of revolu- 
tion are to be considered it is natural to introduce a 
prolate-elliptic coordinate system in the (z, p) plane. 
The equations of transformation from the coordinates 
(z, p) to the prolate-elliptic coordinates (f, 17) are: 



z=2a cosh f cos 7jl 
p=2a sinh f sin jjJ 



(1) 



where O^f and 0<7j<2tt 

Thus f= constant and ^constant represent confocal 
ellipses and hyperbolas, respectively, the distance 
between the foci being 4a. 
For any point in space P(x, y, z) then 



x=2a(X 2 -l)*(l— M 2 )* cos 8 
7/=2a(X 2 — — m 1 )* sin 8 
z=2aXp, 



(2) 



■where \=cosh f and p..=cos jj. 

If, furthermore, the fluid is incompressible the 
velocity potential <j> satisfies Laplace's equation 
A 2 $=0 and since the (X, p., 8) system of coordinates is 
an orthogonal one, takes the form: 



, / 1 , 1 yy n 



(3) 



FLOW PARALLEL TO THE AXIS OF SYMMETRY 

In this case the flow is the same for all meridian 
planes (z, p) and therefore the velocity potential <t> is 
a function only of X and p.. Equation (3) then reduces to 

If this equation is to be satisfied by a product 
*=£(X)M(p.) 

it follows that 

1 dr c£L(xn_ _L_d_Y ( ^ dM(jiT \ 
I(xTSxl (x ~ 1] ~d\ J--M(^L (1 ~^-rfM-J 

which separates into two ordinary differential equations 



d 



(1- 



(5) 



where c is an arbitrary constant. 
Furthermore, if c=n(n+l), each of these equations is 
of the Legendre type and therefore the general solution 
of equation (4) is 



^=^„P„(X)P,(p) 
n-0 



(6) 



This expression for 4> has a singularity at infinity 
since P»(X) is a polynomial of the nth degree in X and 
is therefore infinite for X= » . Since the region outside 
a surface is to be considered and since it must include 
the region at infinity, another solution for L(\) is 
required. This solution, linearly independent of 
P»(X), is the zonal harmonic of the second kind and is 
denoted by Qn(X) where 

It vanishes for X=<» but has a singular point for 
X= ± 1 where it is infinite like log (X± 1). 

Thus, for example, since P 0 (X)=1, P 1 (X) = \, it is 
found that 

n m P" d\ 1, X+l 1.1.1 , 
where | X| >1 and 

n % r- d\ 1 . , x+i 1 
CW=xJ x x^3i) = 2 xl °gx=i- 



,x X 2 (X'-1) 
3X iT 5X 



It may also be shown that 

Q.(X)=|p,(X)log|±f-^(X) 
where jS"„(X) is a polynomial of the (n— l)th degree 
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Another useful expression for <2„(X) is that due to 
F. Neumann (reference 5); namely 

Expanding in decreasing powers of X, 

Expressing X^ in terms of zonal harmonics (refer- 
ence 5) 

i,h 

Xl ~ lj[2i-2le+l] [ikf***™ 
i_ | 

2 I • • 

_ j > used depends on whether % is 



where the upper limit 



(odd and where [27l] = 2 " 4 ' 6 - ' ' 2n; [2n ~ 1] = 1 ' 3 ' 5 ' ' " 
(271-1); [0]=[1]=[-1]=1. Also (2ti-1) /=[2»-l] 
[271-2]. 

Substituting this expression for Xj* in the foregoing 
equation for Q„(\) it follows that 
i u 

• V 2 

^ (x)= i2^ i 2^ 2 ^^^- iP,,(Xi)P< " 2 * (Xi)dXi 

The zonal harmonics P„(Xi) are orthogonal functions 
and satisfy the following relations: 



1 0 if r^s 

P T P4\= 2 if r=8 
2r+l 



Expanding the preceding expression for <2*(X) with 
regard to i and writing the terms with equal indices of 
k in columns and adding these columns, there is ob- 
tained, using the orthogonal property of the P„(Xi)'s, 
the following equation: 

CD 

t-o 

This expression is convergent for |X|>1 and divergent 
for|X|^l. 

Instead of being given by equation (6) the velocity 
potential is now given by the following expression: 



CO 

n-0 



(10) 



which gives the general solution of equation (4) for 
regions outside a surface of revolution and extending 
to infinity. 



In cases of rotational symmetry Where the lines of 
flow are in meridian planes, it is convenient to intro- 
duce Stokes' stream function ty. This function arises 
from the statement that the fluid is incompressible 
(equation of continuity) and is related to the velocity 
potential <j> according to the following equations: 



dip d<p , d\f/ d<p 



(ID 



The lines ty— constant represent the streamlines. It 
may be remarked that, unlike the two-dimensional case 
where both the stream function and the velocity po- 
tential satisfy Laplace's equation, Stokes' stream func- 
tion does not satisfy it. 

The introduction of the variables X, y. into equations 
(11) by means of equations (1) leads to the following 
relations: 

£*0- < fl»t„aJt~ 2„(X.-1)|4 (12, 

If a substitution is made for $ from equation (10) and 
P n (ji) is replaced by its value obtained from Legendre's 
differential equation, that is: 



it is found that 

GO 

, = 2a(l-^(X*-l)2^ f> f +constant (13) 

n-l 

Furthermore, if the body of revolution is moving 
with a velocity U in the direction of the axis of sym- 
metry z, it may be conveniently supposed to be at 
rest and the fluid to have a translation — U super- 
posed on its actual motion. This consideration adds 
a term 2aU\ii to the velocity potential and 2a 2 Z7(l — /* 2 ) 
(X 2 — 1) to the stream function. Therefore 



^=2a 2 (l-M J )(X 2 -l)D 



CO 

l_n-l 



A, dP n dCU 
Un{n+1) dfi d\ 



■] 



+ 1 (14) 



At the surface of the fixed body of revolution the 
normal velocity of the fluid must be zero and therefore 
the boundary must coincide with a streamline 
f = constant, say 0. Hence the boundary condition at 
the surface is given by 



CD 

^J7l 



A n dP n dQ n 
(»+l) dn dX 



-aU=0 



(15) 



In order to find the velocity components U\, u P) in 
the directions of the coordinate lines X, y., respectively, 
it is to be noted that since the system of coordinates is 
an orthogonal one, 



where 



d<j> , d(j> 
ds 1 =dz l +dir=dsx 1 +d8/ 
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By means of the equations of transformation (1), it is 
found that 



1 /X 2 -l Vd* 
^~ 2a\X 2 -W^ 



and 

Therefore: 

and 
Hence: 

SINK-SOURCE DISTRIBUTION 

The distribution of sinks and sources is assumed to 
He along the segment of the axis of symmetry z, 
—2aS z^2a, and to be of intensity _T(zi) per unit 
length. At any point (2, p) in any meridian plane the 
velocity potential due to this distribution is given 
(reference 6, p. 60) by 



4,-- 



' rj-2a[( z — 



"4«-J_ a ,[(a-2,) 5 +p 1 J» 



(17) 



For points lying on the z axis but outside the distribu- 
tion, the velocity potential is given by the simplified 
expression 

4fl-J_2a Z—Zi 

Substituting for z and z u 2a\ and 2a\ u respectively, 
the preceding equation takes the form 



<f> 



1 p X \^,(l/\lJU A, 

4nJ _i X— Xi 



Finally, substituting for <f> from equation (10) and 
noting that P„(l)=l for all values of n, 

n-l 

This is an integral equation for the unknown function 
i"(2aXi). It may be solved in the following manner: 

From F. Neumann's expression for <2n (X) given by 
equation (8) the following development is suggested 
for the distribution function: 

00 

I(2a\ 1 ) =^a»P»(X0 

where — l^Xi^l 

It then follows directly from equation (8) that 



for all values of X. 

Hence On=2vA„ 



and 



J(2aX 1 )=2^VAP»(X 1 ) 



n-l 



(19) 



Thus, given the potential function 4>, that is the 
A n 's, this expression deter min es the equivalent sink- 
source distribution. 

FLOW NORMAL TO THE AXIS OF SYMMETRY 

The differential equation for the velocity potential in 
the case of transverse flow is given by equation (3). 
Recalling that this expression is Laplace's equation in 
the coordinates X, n, 0, it may be solved by supposing 
4> to be a product N{\, n) R(6). Replacing <j> in equa- 
tion (3) by this product, the following pair of differen- 
tial equations is obtained: 



— F- 



X ! -V 



9/*. 
N=0 



(20) 



(21) 



(1-^(X ! -1) J 
The general solution of the first equation is given by 

R=A cos k 6+B sin k 6 
where A and B are arbitrary constants. 
Putting N(\, n)=L(\) M{p) in the second equation 
leads to the following pair of ordinary differential 
equations: 

*[«-*gK*-A>-« 
£[»-<K-£>>«. 

where c is an arbitrary constant. 

Both of the latter equations are of the form of the 
differential equation for the associated Legendre func- 
tions provided, that c = n {n+ 1 ) . A ccordingly , 

M(m)=P„*G0 andi(X) = P„*(X) 

where, for example, 

P.*G*)=(i-/*T ^* 

The general solution of equation (3) may then be 
written as 

^^SS^WP.'IX)^ cos kO+B* sin kO] 

n-0t=0 

This expression, however, has a singularity at infin- 
ity and since only the region outside a given surface of 
revolution is of interest, the infinite region, or the 
neighborhood of X= » , must be considered. Therefore 
P„*(\) is replaced by the associated Legendre function 
of the second kind Q B *(X), where by definition, 



&'(x)=(x*-ir^|^ 



Then 



*i=S cos kd+Bn* sin kO] (22) 

If the body of revolution moves with a uniform 
velocity V in the direction of the x axis, it may be 
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supposed to be at rest and the fluid to have a transla- 
tion — V superposed on its actual motion. Then 



4>=4n+xV 



(23) 



Consider the body profile in any one of the meridian 
planes 0. At any arbitrary point of it the normal 
derivative of <f> is given by 

d<f> j d<p , d<f> j 



(24) 



Since the normal velocity along the meridian curve is 
zero, it follows from equations (23) and (24) that 

(t+^>--(t+^>-° 

Also, x=p cos 0, so that 



dx .dp , 3i a dp 
=cos 6-kt and -3— =cos 0^— 



9«x 



Therefore 



(t+ FMS6 t>'=(t +ycos9 t) fe (25) 

In order that the condition of no flow normal to 
the body of revolution be valid for all values of 0, 
there must be chosen from among all the solutions 
given by equation (22) that one which has cos 0 as a 
factor; namely: 

♦i-S^WWflL'W cos 0 



or 



CO 

0 1 =2aycos0(X 2 -l)Kl-M ! ) 1 2o.^§ 5 (26) 

where 0,=^ 
Furthermore 

ds\=2a Q^2!^ y dX and ds„=2a Q^ ~^Q *<fr* 

so that equation (25) becomes 

^fa+pVcosO) 

Finally, by means of equation (26) and the differ- 
ential equation for the Legendre polynomials the fore- 
going boundary condition takes the following form: 



CO 



an an a\ 



dn 



(27) 



DISTRIBUTION OF DOUBLETS 

The doublets are assumed to have their axes in the 
x direction and to lie along the segment — 2a.gz,g2a 
of the axis of symmetry z. The velocity potential at 



any point (z, p) of some meridian plane 0 then takes 
the form (see reference 6): 

p cos 0 C 2 * J{z{)dzi 



#1 



9 r ta jo 

J_2o[(z— £ 



4r J-2«[(3-2i) 2 +P 2 ] i 
where J(zi) is the intensity of the doublets per unit 
length. 

Substituting for fa from equation (26) it follows that 

v\ di\dQ n _j_ r» j( Zl )d Zl 

V /j * d? 1\ ~^J-^[(z-z 1 ) i +P 1 ]i 

For points lying on the z axis but outside the distri- 
bution this equation takes the following simplified 
form: 



(n+l) a dQ n _ 1 



d\ 16xa 2 



p t 7(2gx 1 )^x t 
J. (x-xo 3 



(28) 



where z x is replaced by 2aX 1} z by 2aX, and (^j~f) 
by n ( n £^\ This is an integral equation for the un- 
known function J(2aX t ). In the solution of this inte- 
gral equation it is necessary that a development of 

-7T — ^ a series of Legendre ploynomials in Xi be 

obtained. The form of this development is suggested 
by Neumann's equation (8). Thus assume that 



C=xT = 2 6s ' Pa(Xl) ^ (x) 



Then substituting this expression for r — r- in- Neu- 

A — Aj 

mann's equation and making use of the orthogonality 
relations satisfied by the Legendre polynomials, it is 
found that 

b H =2n+l 

Therefore 

CO 

^- = 2)(2«+l)Pn(X 1 )$ a (X) 

Differentiating this last expression once with regard to 
X and once with regard to X I} it follows that 

_J = ^n+ldP^dQ^X) 

(X— XO 3 Zj 2 d\ d\ 

Equation (28) then becomes 

n-1 

CO 

It is now obvious that the following assumption must 
ba made: 



J(2aX.) = -8*a*V(l - X, 2 )^,, 

m-l 
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and substituting this expression in equation (29) it 
follows, that 



In order that this equation be valid for arbitrary val- 
ues of X, 

<V=Or 

and therefore 

J(2a\ 1 ) = -8TraW(l-\fi^C a j^ (30) 

Thus, given the velocity potential <j> u that is, the 
O n 's, this expression determines the equivalent doublet 
distribution. 

DETERMINATION OP THE COEFFICIENTS A, AND C, 

Any symmetrical profile may be represented by a 
power series in ji(=cos 17). That is 



X=^o rM ' 



(31) 



r-0 



The rapidity of convergence of this series depends, 
however, on the choice of origin with respect to the 
profile. Since X=ar fl defines an ellipse, the rapidity of 
convergence of the foregoing series may be looked upon 
as a measure of the resemblance of the profile to the 
ellipse X=a 0 . The proper choice of origin may be 
attained in the following manner. The radius of 
curvature R of an ellipse at the end of its major axis 
is given by 

where A and B are its semimajor and sennminor axes' 
respectively. 

Eliminating 17 from equations (1), the following 
equation of a system of confocal ellipses results: 

(2 a cfshrr + (2 a srnhr) 2 = L ^ di8tancebetweei1 
foci is 4a.) 

In terms of elliptic coordinates then 

sinh 2 £ 



B=2a 



cosh £ 



Furthermore, for an elongated ellipse the semimajor 
axis 2a cosh £ is large compared to the semiminor axis 
2a sinh £. This limitation means that £ is small. 
Neglecting powers of £ higher than the second it 
follows that (see reference 7) 

5=2a£ s (approximately) 
The~ends of the ellipse are at 

±2a cosh £=±2o(l ) =± ( 2a +f) 

(approximately) 



and therefore the focus of an elongated ellipse very 
nearly bisects the line joining the end of the semimajor 
axis and the center of curvature. Thus the proper 
choice of origin is the point bisecting the line of length 
4a extending from the point midway between the lead- 
ing edge and the center of curvature of that edge to a 
point midway between the trailing edge and the center 
of curvature of that edge. Having thus chosen a 
reference frame (z, p) in which to present the profile, 
the next step is to obtain the series equation (31). 
This equation may be obtained with the help of the 
following expressions. From equation (1) it can be 
found that 



Hi 



"=2 



[V(r^)"+(&)WGH'+(£)". 



(32) 



where — 1±= ^ 1. 

A series of corresponding values of X and /* are thus 
obtained. In order to express X as a polynomial in ju 
of, say, degree n, it is most convenient to employ the 
method of least squares for determining the (71+ 1) 
constants a r (reference 8). 

FLOW PARALLEL TO THE AXIS OF SYMMETRY 

The boundary condition for this type of flow is given 
by equation (15). In that expression functions of the 

type appear and these are to bo expressed as power 

series in p. 

Suppose the meridian profile to be given by tho fol- 
lowing analytic expression: 



X=Oo+M / fli. 



(33) 



Then on the profile, may be looked upon as a 

function of y. and can be developed in a Taylor series 
in the neighborhood of /x=0 or X=a 0 . That is, 

dQ n _\\ {\-a 0 y d>+ 1 
~Zj Pi 



oV 



p-0 



■■Oo 



Substituting for X— a 0 from equation (33), it follows 
that 

CD / CO \ P 



In the following the expansion of 5" in powers of n is 
to be determined (reference 9, p. 122), p being any 
positive integer and where 



00 

«-0 



Thus 
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or 



5=0 3-0 



CD CD 



r=o p=r 



Expanding S* with, respect to r and writing the terms 
with equal indices of p in columns and adding these 
columns, 



where 



r-0 



In a similar manner, 



where 



where 



p-0 



where 



Hence 



or 



p-0 



V , Oifr^O 

g*,p= / fh,p-r di-i.r and Oo,r= 1 y r=Q 

r-0 



p=0 8=0 



P"0 Q—p 



Expanding according to 2? and writing the terms with 
equal indices of g in columns and adding these columns, 



Q-0 p-0 



(34) 



(35) 



and so on. 

The boundary condition also contains terms of the 
dP 

type where P»(m) is the Legendre polynomial in 11 
of degree n and is given by 

where the upper limit for j is 



n n-1 
2' 2 



pn V\ , ^' [2n-2j-l] ' 

J-0 



n- 



x according as n is 



even 
odd 



n-2 n-1 
2 ' 2 



J=0 



Substituting for ^ i and ^ ! the expressions given by 

equations (34) and (37) into the boundary condition 
equation (15), 



n-2 n-1 



n-1 



>0 



g-0 p-0 



Expanding according to powers of » and writing the terms with equal indices of j in columns and then adding 
these columns, 



.'=0 n-2/+l 



?-0 p-0 
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Putting n—2j—l+g=m this becomes 

vvv... 



\ \ \ 1 „, Uj Asi+i+m-, [2j+l+2m-2g] \^ a p . < ,_ p d p+1 Q ii+l+m . 

ZjZjZj* 1 } (2i+l+m- 2 )(2i+2+m- 2 ) (m-q)![2j\ p! daf+i 

Jz=0 j=0 oj-5 j><=>0 



-aU=0 



Expanding with respect to g and -writing the terms with equal indices of m in columns and adding these columns, 

Vu-YYf IV ■ ^W"-, [2i+l+2m-2 g ] ^ o T ,, g _ p rf^ i+!+m _ 

2-l2-J { 1} (2j+l+m-3(2i+2+m-2) (m-gjl[2j] 2-J p! dZF 1 +aU-0 (38) 



m-0 s-0 j*=Q 



p-0 



If this expression is to be valid everywhere on the boundary surface, it must hold for the entire range of p.. It 
follows that the coefficients of the various powers of p are identically zero. Finally, the introduction of k and 
7i by means of the substitutions q=m—k and p=m—n, respectively, leads to the following expression of the 
boundary condition: 



VV a a -n.n-t 1V [2j+l+2k]A lw+t (fr-'+'Qw+t 

ZjZjk!(m-n)!Zj^ J (2j+l+k)(2j+2+k)[2j\ da*"-™ ~~ 



-aUS 0 m 



(39) 



u-0 i-0 



where 



j=0 



. _ 1 if m=0 , - , . 
5 " 0 if m^O and m=0 » 1 ' 2 > 



Equation (39) represents a set of linear equations 
infinite in number and containing an infinite number 
of unknowns A n . It provides a formal and rigorous 
solution of the problem of potential flow about a 
body of revolution, parallel to its axis of symmetry. 

In the foregoing equations the only unknowns are 
the A n 's. The a*, 5 's are related to the coefficients of 
the power series of X in ju (giving the meridian profile 
equation (33)) and are evaluated by means of equa- 
tions (35). Finally, the Q B 's and their derivatives are 
well-defined Legendre functions. 

For example, if the meridian profile is an ellipse X=Oo, 
then equation (39) becomes 

CD 

W A 3m+m [2j+2m+l] dQsj+i+m 

2J- 1; (2i+l+m)(2i-r-2+m) m![2j] da* 

i ~" =-S a "aTJ 

For m=l, 2, 3, . . . this is an infinite set of linear 
homogeneous equations for the unknowns A3, A3, 

and, since the determinant of the coefficients 

is different from zero, the only solution is that A2, 



A3, are zero. From the first equation, 

i. e., m=0, it is then found that 



Hence 
<*>=- 



. _ 2aU_ 

— ~mr' 

da 0 
2aU 



2aU 



1 w a o+l 

2 log 



|log^ 



Oo— 1 a 0 2 — 1 



a n — 1 Oo*— 1 



If A and B are the semimajor and semiminor axes of the 
meridian ellipse and e its eccentricity, then 



2a— Ae, a„=-; 2a(af— l)l=B 



so that 



UA 



1 1 1+e 



This result agrees with the well-known expression for 
the velocity potential of a prolate ellipsoid of revolu- 
tion (reference 10, p. 132). 



FLOW NORMAL TO THE AXIS OF SYMMETRY 



The case of flow normal to the axis of symmetry will now be treated in a manner similar to the case of parallel 
flow. The boundary condition is given by equation (27): 



From equation (33), 



CO 

^■=a 0 +Y l (n+2)a lin ^ 



(40) 



or 
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Keferring to equation (38) it follows that 

n«l m-0 { n 0 j=0 p-0 

n-»0 m-0 2-0 j— U p-0 

» co tn co q 



n=l ra-0 p=0 



Analogous to equation (34) 



~3^ x (41) 



«=0 p-0 

Hence, in a manner similar to the derivation of equation (38) 

co 

n-1 

« ffl co ft 

^"EE ( - i> , ' : ' ; ' + "7^w m '" R '--' + ^-^ 1 -2/T <«' 

m=i j=o j=o p-o 

Substituting equations (40), (41), and (42) in the boundary condition and equating the coefficients of the various 
powers of ft equal to zero, the following set of equations is obtained: 

CD 

2(- 1 )" [? ^{^H[^^-(2 W +l)(27 l +2)^ +1 ]-27 MIl ^ i .^}=a 0 for h=0 

n=0 

and (after rearranging as for equation (39)) 

A-l k n 



2_i ZjmKk-n)^-"-^-^-"^ - } m ° 2 " 1+ "* da*-»^ 

■0 n=0 m=0 j=0 

4.„VVW"Vf lV [2j+l+2m] ^ 



n-0 m=0 



(43) 



« I n V V a^-n.n- n Yr , v (2j+m) (2j+ 1 + m) [2j- 1 +2m] ,-, ff-'+'tjy,.. > 

n=0 m=0 j=0 

whore A = l, 2, 3, °°. 

This set of equations represents a formal and rigorous solution of the problem of potential flow about a body 
of revolution with flow normal to the axis of symmetry. The only unknown quantities are the infinite number 
of O n 's. The other quantities appearing in the equations are determined as in equations (39). 

If the meridian profile is an ellipse X=at), the <ii,»'s are all zero and, from the second expression of equa- 
tion (43), 

Y [2n+2h+l] \ dQ* +l+h ( 2n-h-2)l 1 _ n 

71=0 

where h—1, 2 ». 
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This expression represents an infinite set of linear 
homogenous equations in the unknowns &, O3, . . . . . 
and, since the determinant of the coefficients is differ- 
ent from zero, it immediately follows that the only 
solution is Gi=Cz=Oi—. . . .=0. From the first 
expression of equation (43), it follows that 



JOT 



•2<3i(ao) 



This result could have been easily obtained from the 
general expression for the velocity potential given by 
equation (26). Thus, assume that 



^!=2ay cos 0(X 2 



-D*(l-MW^f 



If the body moves in the positive direction of the y 
axis with constant velocity V, it may be supposed to 
be at rest and the fluid to have a translation — V super- 
posed on its actual motion. Accordingly 

<j,=2aV cos 0(X s -l)K(l-^w(a^+l) 

At the surface of the ellipsoid of revolution gener- 
ated by the ellipse X=a 0 , the normal velocity of the 
fluid must be zero. Therefore 



-200^+2^(00) 



— Oo 



d<j> 

From Legendre's equation 
Hence (see reference 10, p. 133) 

In the appendix an application of the boundary con- 
dition (equations (39) and (43)) for axial and trans- 
verse flows, respectively, is made to a body of revolu- 
tion obtained from a symmetrical Joukowsky profile. 

INERTIA COEFFICIENTS OF BODIES OF REVOLUTION 

It is of some interest to obtain the coefficients of 
inertia for axial and transverse flows and also to com- 
pare them with those of an ellipsoid of revolution of 
equal fineness ratio (references 11 and 12). 

"When a body moves in a fluid at rest at infinity the 
total kinetic energy of the fluid is given by 



=2aV cos •(^"[ftC^D^+ft W]£-0 



2T=- -af/t^dS 



(44) 



where <j> is the velocity potential of the fluid motion, 
the normal derivative of <j> where the positive direc- 
tion of the normal to the surface of the body is into the 
fluid and the integration is performed over the surface 
of the body; <r denotes the density of the fluid. 



FLOW PARALLEL TO THE AXIS OF SYMMETRY 

Since the velocity potential of this type of flow is 
independent of the angular coordinate 0, the following 
equation may be written for the element of surface: 

dS=2T P ds 

where ds denotes the element of length along a merid- 
ian profile. Hence, 



2T= —2Tra<p P 4> |^ ds 



If the body moves in the direction of its axis of sym- 
metry with a uniform velocity U the boundary 
condition is 

^-ds— — Ujrds 
on on 

Also, according to equation (24) 
Therefore, 

2T= -87raV uj' <f> [(1-/X 8 ) X^- (X 2 - 1)m](Zm 
In general then 



2T= -S^roVC/Jl [(1 -m 2 ) X ^- (X s - l)M^A n PM Q n (\)d» 



(46) 



If M is the mass of fluid displaced by the body, then 
the coefficient of inertia k a is the quantity multiplying 
MZP in the expression for 2T. 

If the body is a prolate spheroid X=Oo the foregoing 
expression for 2T becomes: 



2T=|™(2a) 3 £7 2 (a () 2 -l) 



1 , gp+1 



l 1no ,gQ+l 

a 0 2 -l 2 10g ao-l 



But 2a=Ae, ao=— and 2a(a<?— 1)*=B where A, B 

6 

are the semimajor and somiminor axes, respectively, 
and e the eccentricity of the elliptical meridian section. 
Therefore 

1 , l+e 1 
2T=—=£ r %TrAB i <rU 1 =k a MU 2 
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The coefficient of inertia for a prolate ellipsoid in 
uniform axial motion is then given by 



1+e 



k a = 



2e 



1 1 , l+« 



(47) 

(See reference 10, p. 144) 

Equation (46) is now evaluated, for the case of a 
body of revolution obtained from the Joukowsky pro- 
file 61=0.15, €j=0.10. (See appendix.) The volume 
of this body is found to be 

Q=|ir(2a) 3 X0.05342 

so that the expression for 2T may be written: 
^ 2,(2^3X0-003139 ^^ 
|*-(2a) 3 X0.05342 

or jfc a =0.0881. (See table I.) 

Compare this value of k a with that of a prolate sphe- 
roid whose fineness ratio is the same as for the above- 
mentioned body of revolution. The fineness ratio / is 
denned as the ratio of the length to the maximum 
diameter of the body. The maximum diameter is 
obtained from equation (53) by means of the condi- 
tion ^=0 and the length of the body is given by 

Z=2a(X,, Dl +Xp__i). By means of these expressions it is 
found that/=4.208. The fineness ratio for an ellipse 
is given by 

, A 1 



or 



■V 



1-^=0.971 



where e is the eccentricity of the ellipse. 

Substituting this value of e into equation (47), the 
following value of k a is obtained, 

fc a =0.0757 

A theorem enunciated by Munk (reference 13) states 
that when the disturbance caused by a body moving in 
an infinite fluid is replaced by fictitious sinks and 
sources, the total mass is the sum of the products 
obtained by multiplying the intensity of each source or 
sink by the potential of the parallel flow. This 
theorem will now be shown to be only a first approxi- 
mation and to hold exactly only for ellipsoids of 
revolution. Thus from equation (19), 



n-l 



where Zi=2aX, 



The strength per length dz\ is then 

CO 

iirdffS^A* P»(Xi) d\ 
n-l 

The velocity potential at (z u 0) of the parallel flow is 
given by 

<j>=2a 

Hence, according to Munk's theorem 

2T tolal =8vaa i U^ i A H j\ 1 P,(X0 dh 



or 



Therefore 



2T tolal =^(2aycrZP^ 



2 r^-(|^ 2 4-i)w 

where Q is the volume of the body and M is the mass 
of the displaced fluid. 

The coefficient of inertia for axial flow is therefore 

. _±ir(2a) 3 A 1 

Ka ~3~^~2aU 

This expression lor k a is valid for a prolate ellipsoid 
but is not valid for a more general shape. 

It is obvious from this expression that Munk's 
theorem applies exactly only to ellipsoids of revolu- 
tion since only the coefficient A x appears. 1 In order to 
provide a numerical comparison between Munk's 
theorem and the exact method, the foregoing equationis 
evaluated for the body of revolution whose meridian 
curve is the Joukowsky profile ej=0.15, es=0.10. It 
yields a value & o =0.0717 as compared with the more 
exact value k a = 0.0867 already obtained by means of 
the fundamental equation (46). 

FLOW NORMAL TO THE AXIS OF SYMMETRY 

For flow normal to the axis of symmetry the velocity 
potential depends not only on the elliptic coordinates 
X, p, but also on the cylindrical coordinate 6. Hence, 
the equation for the element of surface dS is 



dS=pd6ds 



and equation (44) becomes 



2T=-<rSS P 4>fadsd3 

If the body moves in the direction of the transverse 
axis Ox with a constant velocity V the boundary 
condition is 



= — V cos B^ds 



1 The above criticism of Munk's theorem has been found to be incorrect. (It Is to be noted that the volume of the body contains all of the coefficients At Implicitly.) 
This theorem may be readily verified by applying Green's second theorem to the space Internal to the shape and enclosing the appropriate distribution of sinks and 
sources-)'. Then If*]- Kr and (where Aj*+7—0),Qreen'ssecond theorem together with equation (44) Immediately yields the following expression on Munk's theorem: 



2T/i.u 



+ TP J> (volume of body) ■= 
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Also, in general 



CO 

*=2* V cos 8Q?-1)* H-^O/^ 



Hence 



For a prolate spheroid, X=a 0 and 

4,=2aV cos 0(X 2 -1)J (1-^*)* 

where 



From equation (45), it follows that 

^ds=-8a 2 (l-/x i )i (X s - 1)1 (jxd\+\dix) 



0,= 



OS 

2r=(2 a) 3^£(l_^(X»-l)(,|+x)2af"^^ (48) 

n-l 

00 

J( 2l ) = -8^ 2 F(1 - X^JJ C, § 

n-l 

Then according to Munk 

2T tota i=<rVj^ e 7(2aX,)2oiX 1 



ig-2<M«o) 



Therefore 



2T= 



ao 2 -! 2 10g Oo-l 4 



goi.- g<>+ 1 , ^ o 3 

2 log 5^I + 57=T- 2 



or 



1 



(49) 



6 s ^ 2«» 



>l-e 



(See reference 10, p. 145.) 

For the body of revolution whose meridian curve is 
the Joukowsky profile e 1 =0.15, ^=0.10 (see appendix) 
it is seen from equation (48) that 



2T= 



7 ro-y i (2a) 3 X0.059587 
| ir(2a) 3 X0.05342 



Q=0.8366XMF 2 



Therefore k T = 0.8366. (See table II.) 

According to equation (49) for the prolate ellipsoid 
of equal fineness ratio/=4.208 and £ r =0.8689. 

According to Munk's theorem the inertia coefficient 
k T of transverse flow may be obtained from the doublet 
distribution along the axis of symmetry. Again, as in 
the case of axial flow, this theorem is a first approxima- 
tion and holds exactly only for ellipsoids of revolution 
since an expression for k T is obtained that contains 
only the coefficient &. Thus from equation (30) it 
follows that 



or 



smce 



to 



=1 



Hence 



and 



or 



3x7 



2T Mal =-%™{2 a fV i C x 



| ff (2a) 3 C l +« 
2W=- ^ MP 



8 



3 

In order to give a numerical comparison between 
Munk's theorem and the exact method, the fore- 
going equation is evaluated for the body of revolution 
whose meridian curve is the symmetrical Joukowsky 
profile ei=0.15, 62=0.10. It yields a value Zr r =0.8210 
as compared with the more exact value A>=0.8366 
obtained from the fundamental equation (48). 
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APPENDIX 



APPLICATION OF THE ANALYSIS TO SURFACES OF 
REVOLUTION OBTAINED FROM SYMMETRICAL JOU- 
KOWSKY PROFILES 



By means of the mapping function 



(50) 



the circle ki of radius a in the f plane is transformed 
into the line segment (—2a, 0; 2a, 0) in the f plane 
and the circle k 3 of radius (l+ei+e^a with center at 
(ei a, 0) is transformed into a symmetrical Joukowsky 
profile J in the f plane. (See fig. 2.) 




e.a 



Fioube 2.— Mapping of a circle Into a symmetrical Joukowsky profile. 

If in the f plane PQ=ae^, PO=aei, angle POQ=v, 
and angle PQs'=<f> then, according to the law of 
cosines, 

gS«-W = i-j-25coscS+5 ! (51) 



where 



S=- 



(52) 



"1 + 6! + ^ 

Again, by the law of sines 

sin <j> 
tan v = i_i_„ n „ j. 
5+cos <f> 

Putting f / =oe t+ft into equation (48), 
f=2a cosh (£+iy) 

or z=2a cosh f cos tj, p—2a sinh £ sin 17 

The latter two equations are, in fact, the equations of 



transformation from the rectangular coordinates (2, p) 
to the elliptic coordinates (£, 77). Since (2, p) refer to 
points of the Joukowsky profile J, using equations (51) 
and (52), the following parametric equations of the 
system of symmetrical Joukowsky profiles may be 
obtained 



,_ 6 1 (i+2g y +a 2 )^ 



25 
S+v 



2£ 1 (l+25f+5 2 )» 



'(1 + 25H-S 1 )* 



(53) 



where X=cosh £, p.=cos 17, and v=cos (the inde- 
pendent parameter). 

From these equations X can be expressed as a power 
series in n by means of a Maclaurin's expansion in the 
neighborhood of m=0 (i. e. v=— 8). Thus, 1 



X=o 0 (l + a T +^-^+l 7 B...) 



(54) 



where 



6l 2 (l— 5 2 )+5 g g 2 )— 5 a , 

ao= 2e 1 5(l-« 4 )«' a= V(l_5*)+5* 011(1 



FLOW PARALLEL TO THE AXIS OF SYMMETRY 

Equation (39) is a set of linear equations for the 
mfinite number of unknown coefficients A n and pro- 
vides a solution of the problem of axial potential flow 
about a body of revolution. In practice it is necessary 
to evaluate only the first few coefficients A n . From 
equation (39), neglecting the A n 's after A s , the follow- 
ing equations are obtained: 

' This power series suggests the form 

for n from eqoaUon (53) Is herein substituted the equation 



In fact, If the 
forX Is obtained 
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a **$&\ A -o 



+2 



Jj>(5a 1>0 — a M ) 



"Sao* 



8[24a M ^+12(2 ai . 0 a 1 ^+ ai , i )^+12 ai , ai ^+ ai ^]^ " 
+32(6^ +6a,oa u ^+a 1 ^> a 

+2[24(5a ia -a 1 ,,) ^ +12(5^^-2^^^-^^^- 12^,^.0^-^^^]^ 

+ 24[2(7a 1 ,-3a u )^-6a 1 . 0 a, 1 g -^^]^ 

+ [5(M^-24(14a u -a I ^+12(2a 1 .oa 1 ^+a u s -14a 1 /)^ + 12a 1 ,o 1 /^+a 1 .o^]A=0 



(55) 



The coefficients of the unknown A n 's can be calcu- 
lated, simply by a knowledge of power-series develop- 
ment of X in (i- e., the quantities do/Oi^, a u , . . . are 
obtained from equation (33)). The zonal harmonics 
Qn (oo) are given by means of the recursion formula: 

(n+ 1) Qn +1 (oo) - (2n+ 1)^^) +n&-i (a*) =0 



In the apphcation of this recursion formula it is neces- 
sary to calculate Qo(oo) and & (<*<>) independently where 

Qo(oo)= |log g±j and &(ao)=0o&(a<,)-l 

The first derivatives are then obtained from the 
following relation: 



^+g = (l+ao)^ + Q,(a 0 )+3(2 I (ao) + 5Q a (a 0 )+ . . . + (2 W +l)&(ao) 



where it is necessary to determine independently. 

In order to calculate the higher derivatives the pre- 
ceding recursion formula may be repeatedly differen- 
tiated with regard to Oq. The higher derivatives of 
<2o and Qi are obtained independently by means of the 
(r— l)th derivative of Legendre's differential equation: 



(oo 2 



daf' 



If the constants Oo, Oi.i> «uj • • • and the various 
derivatives of Q t are known, the coefficients of the 
unknowns A t in equation (54) are easily calculated. 
The resulting system of linear equations can then be 
solved for the A t 'B which in turn determine the poten- 



tial function <j> given by equation (10). A knowledge 
of the potential function yields directly the velocity u 
given by equation (16). Finally, according to Ber- 
noulli's equation, the pressure p on the surface is 
given by 

or 



n __ 7>—Va 



=l-(^where Pj =|Z7 s 



The numerical work is straightforward but somewhat 
tedious owing largely to a lack of tables of the zonal 
harmonics of the second kind. 

As an illustration of the procedure here outlined, 
consider the body of revolution whose meridian curve is 
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the Joukowsky profile (reference 14) defined by ex =0.1 5 
and 62=0.10. (See fig. 2.) From equation (53) there 
can be written then for the power-series development 
of X: 

X=1.02340+0.02630m+0.007476m 2 

-0.0000273^+ • • • , (56) 

a very rapidly convergent series. 



Figure 3 shows a graph of this function with 
Z(2aXi)/4iraZ7as ordinate and Xi as abscissa. 




Figure 3— Sink-source and doublet distributions. 



In order to obtain the pressure distribution, the 
following expressions can be evaluated at a sufficient 
number of points of the boundary: 

©' 

Note here that the 



is exact for the body of revolution obtained by super- 
posing a uniform velocity U on the flow from the sink- 
source distribution given by equation (58). This body 
is a very good approximation to the actual body 
obtained by revolving the Joukowsky profile about 
the axis of symmetry, so that in calculating the pres- 
sure distribution it is permissible to use the (X, /i) 
values as given by equation (55). 

Table IV shows the sequence of operations to be fol- 
lowed in obtaining the pressure distribution and figure 
4 presents graphically the pressure distribution. 



Here 

Oo=l-02340, di o=0.0263, <h !=0.007476, 
«w=0, <hs=— 0.0000273, 

The zonal harmonics of the second kind and their 
derivatives are given by table III. 

Substituting these numerical values into equations 
(55) the following set of five linear equations is ob- 
tained for Ai, At, As, Ai, and A&: 



(57) 



(58) 



5 




n-1 

5 



n-1 

and then substituted into equation (16) for u 2 where u 
is now the velocity at the surface of the body with the 
body considered to be at rest with regard to the fluid. 
The velocity u is calculated by means of the following 
expression: 



/.Or 




-5 L 

i iuuek 4.— Theoretical pressure distribution (axial flow). 



3.780 At +1.400 A s =2aU 

2.342 A^+ 4.333 A+ A=0 
0.810 A 3 - 0.850 At+ 
0.763 A 3 + 0.300 A*+ 



A=0 
A=0 



19.385 Ai 

10.372 ^- 7.635 4 2 - 

- 1.152 At+ 1.294 A t - 

- 0.500 Ai+ 0.685 A t - 

-27.957 4,-1-39.785 4 s -57.432 4,+63.649 A,+ 
The solution is given by 2 

^=0.0573 X2aU, A i =0.0726X2aU, 4 3 =0.0296X2aZ7, A=0.0065X2aZ7, ^=0.0015X2(1/7 
The sink-source distribution obtained from equation (19) is: 

J(2aX 1 )=4ira7J(-0.0339+0.0157X 1 -l-0.0846X 1 i +0.0606X 1 3 +0.0283X 1 4 -f-0.0119X 1 s ) 



velocity potential 

«S=2aC7I0.0573 PMQM +0.0726 P 2 (m)&G*) + • • • +0.0015 Psto&Qx)] 



» In its exact form the system contains an infinite number of equations with an infinite number of unknowns At, At, ... . For practical purposes, however, the 
following method of solution Is suggested. Suppose the system of equations to have been solved to an arbitrary degree of approximation, say three. Then to this solu- 
tion there corresponds a definite sink-source (or doublet, as the case may be) distribution from which can be obtained the corresponding profile and hence a (X, p) 
curve. This (X, /i) curve can then bo compared to the (X, p) curve of the actual profile. In order to Improve the approximation, the true (X, p) curve can be shifted 
In such a manner that a repetition of the process of solution, to the same degree of approximation, yields a new system of (X, p) values closer to the actual set of (X, p) 
values than the first approximation. In this manner the process can be carried on until the desired degree of accuraoy Is obtained. 
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FLOW NORMAL TO THE AXIS OF SYMMETRY 

From equations (43), for the first five coefficients 
Oi, Ci, Cs, Ct, G s the following set of linear equations 
can be obtained: 



A t 1 C l +A t *C2+ AfQ+ASCt+AfC^i a u ,_ a (60) 

where i=2, 3, 4, 5 

and 



A'=-<l<".fH«*) 

Al n (n ^ 9^ 

^ 2 =3[a,o 2 g+( ao +2a l . 1 ) 6ft] 

^- 3 [¥»+«"««"+««»^ +(3 «'--- 7<, ">3&] 

4i , = _ 3 |-»ig + l (aA ,_ 16o ,, l!) ^_l (5at+ 33« 1 .0g + 30«,] 

^-S 5 ¥^+<««<'-.-^ i >0-<"°'+ s "'"'S+ 42oa ] 

^ 2 = 3 [-3f ^+-,o s ^^^+(^.»a 1 ,+2a, 1 2 +^, 1 -4 ai ^ g+<4*...-^og] 
^3 = _^g + | aio(aoai>1 _ 7aiioi ) g+| (a 0 a, 2 -43a 1 ^, 1 -5a 0 a 1 . 0 ) ^+3(35a 1 ..-22a 1 . a) g] 

^=-f[^^+3a,/^ 
^ = |^!g+fa, 0 (a^ A -l^ 

^=^+^-^W 
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+ (501,301,0+ 5ai i2 ai,i+a 1) a 1 o—10o li oa u )^^ + (5o M — 11©!^) 

A 3 =-|ff^+^ (3o,Oi,-13oi^ ^+|[ao(oi, s +2ai.oO,^)-o,., s (5ao+53a u )]^ 

+3[ao(a w -5o u )+2a w (20a 1J -27o 1 ^)-27o ia s ]^ +15(17a u -lla u ) ^ j 
5f" 3a 1J o B ff& . jt 5o u +a 0 d t Q < ^ Oi.qOij+0,.1 8 , 22 
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+ (15a 1J ai, 0 + 15a u a u +3aoOu— 135^^)01,1— 7aoOi^) + (15ai, 4 — 138a u + 161a w ) 



^- 5 (^^+^( 6 ^-71o 1 /)g+^[o 0 ( ffll /+2a 1 ^ u )-Oi/(14a 0 +143oi J )^ 



OOo 3 



+3[oo ai - 8 ~ 14a '- 1 -9(o., 1 »+20i J a u ) + a./]f§+§ (21a 0 -145a„+658a I ,) ft j 

Again, substituting the numerical values given by equation (53) and table III into the foregoing expres- 
sions, there result the following equations: 



-22.402 Oi— 1.363 C 2 + 34.018 0+ 2.578 0<- 39.602 0 5 = 1.02340 
25.036 0-67.392 O s - 45.658 0+ 163.240 0,+ 66.762 O 5 =0.05260 
-20.665 Ci+78.034 O-132.860 0>— 221.940 O+497.330 C 6 =0.02243 
-15.235 O+64.050 0,-178.260 O+202.310 O.+391.330 O 5 =0 
10.437 0-46.495 O+150.380 O-361.150 0+197.350 O 5 =0.00014 



(61) 



The solution is given by 



0= -0.0486 
O=-0.0178 
O=-0.0027 



0= -0.00028 
G s = -0.00005 

The doublet distribution function J(2a\{) then be- 
comes: 



J(2aXi)=87ro J F(l-Xi sl )(0.0447+0.0512Xi+0.0187Xi 2 +0.0049X 1 3 +0.0022Xi 4 ) 



(62) 



The graph of this function with JlSvaW as ordinate 
and Xi as abscissa is shown in figure 3. 

DETERMINATION OF THE TRANSVERSE-FORCE DISTRIBUTION 

When the axial flow is combined with the transverse 
flow some information regarding the distribution of 
forces over the surface of the body can be obtained by 
introducing the notion of the transverse-force coeffi- 
cient. For the pressure difference at the surface, 
according to Bernoulli's equation: 

p-Po=%OP+V>-<?) 

where .p 0 is the pressure at an infinite distance from 
the body and q is the velocity of the fluid at the 



surface, supposing the body to be at rest and the 

fluid to strike it at an angle a where tan a=-y 

. Now g has three components — in the directions dsx, 
dsp, and pdd. They may be denoted by g>, and qe, 
respectively. Also denote by u\ and ii p the velocity 
components of the axial flow and by £x, r p the velocity 
components of the transverse flow taken in the plane 
0=0. Then, 

2x=«x+tfx cos 6 
2p=«p+»d cos 9 
ge=ve sin 0 



Introducing these -values into Bernoulli's equation, 
p— Po=^U 3 +V 2 — u\ 2 — uf— vx* cos* 0— V cos* 0— v$ 2 sin 2 0— 2(uxV\+u ll v f >) cos flj 
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It is only the term 

ffCttxfx+Wpfp) cos d=<r(uv) cos 0 

that need, be considered, for the other terms have equal 
values for 8 and ir-6 and accordingly vanish when 
integrated over the entire cross section. The resulting 
transverse force, relative to an annular element of 
width unity, is therefore 



Elemperer defines the transverse-force coefficient by 
dQ 

(*«0 



0= 



dz 



•=2ttp 



U'+V 2 



or 



(63) 



By means of the velocity potentials of the axial and 
transverse flows this last expression takes the following 
form: 3 



^Q='(rJ\u v) cos 2 6 p de=aTr p (« v) 

L\n-1 / ln-1 

-(22^f t «-w + 0iS a S[ w(w+i ^-° i >-' i f ! ] + ' i ll 

\n-l / U=l JJ 



2av sin 2a 



(64) 



Tables IV and V give the numerical data for the 
evaluation of the right-hand side of equation (61) and 
figure 5 represents graphically these numerical results 



./5r 




Fioube 5.— Transverse-force distribution. 



—. — k~ as ordinate and ~— (= X/i) as abscissa. 



With s • K " ao uiiuuuto <uiu 0 

2air sin 2a 2a 
According to theory the positive and negative areas 
included by the jS curve and the x axis are equal; that 
is, there is no resultant lift force but only a simple 
couple. 
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<fX 
if 




(1-M>)X 
dX 
"A 


-(1-a')X 
dX 

3? 


5 

n-l 


r 

(x'-Dp 

L 

6 

7l n l 


1 


a 11767 


a 04114 


0 


0 


a 11767 


a 10026 


a 01179 


.9 


.09814 


.03968 


.20000 


.00794 


.09120 


.08114 


.00740 


.8 


.03087 


.03821 


.37772 


.01443 


.06624 


.06337 


.00420 


.7 


.06611 


.03673 


.63319 


.01968 


.04653 


.04696 


.00214 


.8 


.06129 


.03626 


.66680 


.02361 


.02778 


.03191 


.00089 


.6 


.03916 


.03376 
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.02629 


.01287 
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.00023 


.4 
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.02824 


-. 02296 
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0 
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.00081 


-.1 


-.00421 


.02481 


L 01064 


.02607 


-.02928 


-.03660 
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-.04132 


.00107 
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-.01217 
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.61441 


.00816 
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-. 04433 


.oooso 
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-.01146 


.01439 


.36267 


.00622 
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-.04203 


.00070 
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-. 01041 


.01292 


. 19110 


.00247 


-. 01288 


-. 03895 


.00050 


-1 


-.00912 


.01146 


0 


0 


-.00912 


-.03542 


.00032 



According to the trapeioldal rule 

5 

Jij [(X»-%-(l- P i)X^]S^^P.0i) Q.(X)a>-0.003139 



n 


Al 1\ 

(L— 11*) 




(X»-l)(l-p>) 

l"d7 fx ) 


5 

n=l 


ftl— 1)(1— B J) 

( p ii +x ) 

J, dP.dQ. 
^ all OA 


1 


0 


L09839 


0 


a 66858 


0 


.9 


.02072 


1.03882 


.02256 


.76363 


.01712 


.8 


.03630 
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.7 
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.6 
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.6 
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.80338 
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.4 


.06004 
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.82192 


.05172 


.3 
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.06168 


.83623 


.05149 


. 2 


. 05641 


L 03482 


. 05837 


. 86079 


.01966 


.1 


.05237 


L02889 


.06388 


.86690 


.04666 


0 


.04735 


L 02340 
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.04363 


-.1 
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L 01837 
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.89411 


.03799 


-.2 


.03673 


L 01378 


.03622 
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.03315 
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.02968 
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.02997 
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.02796 


-.4 


.02382 


L00595 


.02396 


.96040 


.02278 


-.6 


.01828 


L 00269 
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.96979 


.01777 


-.6 


.01327 


.99990 


.01327 


.98811 


.01311 


-.7 


.00886 


.99754 


.00384 


L0OS04 


.00391 


-.8 


.00616 


.99562 


.00613 


L03034 


.00529 


-.9 


.00219 


.99414 


.00218 


L05385 


.00230 


-1 


0 


.99309 


0 


1. 27840 


0 



According to the trapezoidal rale 
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TABLE ILL— (oo= 1.0234) 



n 


<3. 


da% 


dm? 


da? 


d'Q. 

dan* 


<!'<?. 

dad 1 


0 
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-2L1202 


913.16 
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I 
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-19.3846 
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2 
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3 
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4 
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5 
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Pi 


Pi 


Pi 


Pi 


Pi 
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1 
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1 


1 


1 


1 


1 
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0.2072 
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.40 


-.2800 
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TABLE IV— Continued 
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8.646 


-. 0102 


.3 
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22.801 
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0 


—3 


0 
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0 
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31.801 


-.3 


-2.46 


4.3425 
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37.607 
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5.16 
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-L60 
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